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Abstract. In this paper we consider two aspects of the stability of a compressible mixing layer: 
Absolute/Convective instability and the convective Mach number. We show that, for Mach numbers less 
than one, the compressible mixing layer is convectively unstable unless there is an appreciable amount of 
backflow. We also present a rigorous derivation of a convective Mach number based on linear stability 
theory for the flow of a multi-species gas in a mixing layer. Our result is compared with the heuristic 
definitions of others and to selected experimental results. 
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1. Introduction. In the past few years there has been an increasing interest in the sta- 
bility of compressible mixing layers. This results, in part, from the projected use of the 
scramjet engine for the propulsion of hypersonic aircraft. Knowledge of the stability of 
the mixing layers in the combustor may allow one to control the downstream evolution of 
the flow. 

Both temporal and spatial stability calculations have recently been carried out for 
compressible mixing layers. If the flow is absolutely unstable its stability is properly 
described by temporal theory, but if it is convectively unstable spatial theory is appropri- 
ate. It is known that the incompressible mixing layer is generally convectively unstable 
but that a sufficiently large backflow can make it absolutely unstable (Huerre and 
Monkewitz, 1985). Such information is not available for the compressible mixing layer 
for which temperature effects are important. 

The first part of the work reported here is to determine whether a compressible mix- 
ing layer is absolutely or convectively unstable. The calculations are presented in section 
2 and results are given as a function of the Mach number of the fast stream and the velo- 
city and temperature ratios of the streams. 

An additional reason for the renewed interest in the stability of compressible mixing 
layers is that a number of experiments (e.g.. Brown and Roshko, 1974; Chinzei, Masuya, 
Komuro, Murakami, and Kudou, 1986; and Papamoschou and Roshko, 1986, 1988) and 
numerical simulations (e.g., Guirguis, 1988; Lele, 1989; Sandham and Reynolds, 1989; 
and Mukunda, Sekar, Carpenter, Drummond, and Kumar, 1989) have shown that the mix- 
ing layer becomes less unstable with increasing Mach number. In order to correlate the 
experimental results, a number of experimentalists have used a heuristically defined "con- 
vective Mach number". This idea, first introduced by Bogdanoff (1983) for compressible 
flows, has resulted in several different definitions of a convective Mach number. In addi- 
tion, when a mixing layer is formed between two different gases some definitions will 
yield different values for the convective Mach number in the two gases. It is not obvious 
which value is the proper one to use in the correlation, or whether some average of the 
two is appropriate. 

We have recently earned out extensive spatial stability calculations for a compressi- 
ble mixing layer (Jackson and Grosch, 1989a,b). Our results suggested a way to 
rigorously derive from linear stability theory a single convective Mach number for a 
compressible mixing layer for both a single species gas and a multi-species gas. In sec- 
tion 3 we present that derivation and discuss its application. We also present a com- 
parison with selected experimental results. 

2. Absolute/Convective Instabilities. The concept of absolute and convective insta- 
bilities was introduced by Briggs (1964) in the context of plasma instability. The same 
idea was put forward independently by Gaster (1968, 1981) within the context of the clas- 
sical theory of hydrodynamic stability based on the Orr-Sommerfeld equation. Very 
recently these concepts have been applied to classify the instabilities of incompressible 
flows; for example, the mixing layer (Huerre and Monkewitz, 1985); capillary instabilities 
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(Leib and Goldstein, 1986a,b); and wakes (Hultgren and Aggarwal, 1987; Monkewitz, 
1988). A review of these concepts can be found in Huerre (1987). 

A flow is said to be absolutely unstable if the response to an impulse in space and 
time is unbounded everywhere in space for large time. On the other hand, if the response 
decays to zero everywhere in space for large enough time, the flow is said to be convec- 
tively unstable. In this latter case the response to the impulse is a wave packet propagat- 
ing downstream from the source with the waves forming the packet having growing ampli- 
tudes. 

Here we examine the stability of the compressible mixing layer in order to determine 
whether it is convectively or absolutely unstable. We have shown (Jackson and Grosch, 
1989b) that the qualitative stability characteristics of a compressible mixing layer are 
independent of the detailed shape of the velocity profile. Therefore in this study we 
choose to approximate the mean velocity by a hyperbolic tangent profile and to calculate 
the temperature profile from the Crocco relation. Thus 

U [l + P(/+d-Pt/)M^(Tl)]. (2-1) 


7 = 1 - 


1 - Pr 

l -Pw 


(l-U) + l-Y-M 2 (l-U)(U -Vu). 


( 2 . 2 ) 


All quantities have been nondimensionalized by the values of the variables in the fast 
stream, taken here to be at +<*>. Here P^ is the ratio of the speed of the slow stream to 
that of the fast stream, Py is the ratio of the temperature in the slow stream to that of the 
fast stream, M is the Mach number of the fast stream, y is the ratio of specific heats of the 
gas, and T) is the similarity variable in terms of the downstream distance and the 
Howarth-Dorodnitzyn variable. Note that 

— 1 < P u < 1, and py > 0. (2.3a,b) 


It is well known that the instabilities of compressible free shear flows are predom- 
inantly inviscid, except at very low Reynolds numbers. Consequently, we considered only 
the case of inviscid instability. The formulation and the method of solution of the linear- 
ized instability problem for a compressible mixing layer was given by Jackson and Grosch 
(1989a). The eigenvalue is a zero of the characteristic equation relating the wavenumber 
a and the frequency to at fixed Mach number. Since a(oo) has a square root branch point 
singularity at a zero of the complex group velocity 3 to / 3a, transition from convective to 
absolute instability occurs when the zero lies on the real to axis. We therefore choose to 
to be real, a to be complex, and carried out a numerical search for a zero of 3co / 3a. In 
these calculations we took y = 1.4 and varied p^, Py, and M. 

Figure 1 shows results for zero Mach number. The locus of the branch point position 
separating the regions of absolute and convective instability in the p;y - Py plane is 
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shown. In the region to the left of this curve, the flow is absolutely unstable, otherwise 
the flow is convectively unstable. For pj = 1, (3^ = -0.136 which is the value found by 
Huerre and Monkewitz (1985). Note that increasing P 7 - requires a substantial decrease in 

to make the flow absolutely unstable. On the other hand, for a fixed negative {3^, 
cooling the gas of the slow stream can cause the flow to change from convectively to 
absolutely unstable. When (3j is small, the equations become stiff and for this reason we 
were not able to carry out calculations for P r <0.15. Figure 2 shows related results. The 
curves separating the regions of absolute and convective instability in the P u - M plane is 
shown for three different values of P 7 -. In the region below each curve, the flow is con- 
vectively unstable and above the curve it is absolutely unstable. Consistent with the 
results of Figure 1, increasing P^ causes a shift in the curve to more negative values of 
P u at any Mach number. A consistent trend is apparent for all values of P 7 -; as the Mach 
number increases the amount of backflow necessary to cause absolute instability increases. 

The results presented in the above figures shows that, for the standard co-flowing 
mixing layer, the flow is convectively unstable. This has been shown only for subsonic 
disturbances. Because supersonic disturbances have a very different behavior at either +«> 
or - 00 , substantial modifications to the Green’s function is required. Flow control of the 
convectively unstable compressible mixing layer can now be achieved by forcing at the 
frequency of the most unstable mode obtained from spatial stability theory. However, 
Morkovin (1988) notes that the criterion described above for absolute/convective instabili- 
ties is strictly valid only for small disturbance linearized parallel flows. Nonparallel and 
nonlinear effects can alter these results. 

During the course of this study, we became aware of the work of Pavithran and 
Redekopp (1989). In their work they also examined the absolute/convective instability of 
a compressible mixing layer. The velocity profile was given by a hyperbolic tangent, as in 
this study, but the temperature profile is not given by the Crocco relation. Instead, their 
temperature equation is equivalent to that obtained from (2.2) by dropping the Mach 
squared term. Despite the fact that this is not strictly correct, it yields a good numerical 
approximation over the range of Mach numbers considered. Our results differ from those 
of Pavithran and Redekopp by at most 10%. 

3. The Convective Mach Number. The motivation for this study is implicitly con- 
tained in the results obtained from our study of the stability of the compressible mixing 
layer (Jackson and Grosch, 1989b). In this study we took (3^ = 0, used various values of 
(3j, varied the Mach number from 0 to 7, and used a number of thermodynamic models 
with different values of the Prandtl number and different viscosity-temperature relations. 
In all cases, we found that the maximum growth rates of the disturbances exhibited a simi- 
lar dependence on Mach number. As the Mach number was increased from zero the max- 
imum growth rate decreased by a factor of 5 to 10 up to M * , the Mach number at which 
the sonic speeds of the two streams are equal. For Mach numbers greater than M* the 
maximum growth rates level off. This suggests that M* is the proper Mach number scal- 
ing for this stability problem. Thus, a convective Mach number can be defined as 
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M c = M / M* . (3.1) 

With this definition, growth rates decrease for M c < 1, while they level off for M c > 1. 

Motivated by this observation, we derive M* for a multi-species gas as a function of 
the density, temperature, and velocity ratios of the streams and of the thermodynamics. In 
section 3.1 we formulate the problem, and in section 3.2 the results are given. 

3.1. Formulation of the Problem. The problem considered here is that of the invis- 
cid spatial stability of a steady two dimensional compressible mixing layer which lies 
between two streams of N multi-species gases with different speeds and temperatures. We 
first consider the case of two dimensional disturbances, and then generalize our results to 
that of three dimensional disturbances. The nondimensional inviscid equations governing 
this flow are (Williams, 1985) 


p, + ( p u ) x + ( p V ) y = 0, 

(3.2a) 

P ( u t + U u x + V U y ) + C ! P x = 0, 

(3.2b) 

P ( V, + U v x + V V y ) + C ! Py = 0, 

(3.2c) 


p ( H t + u H x + v H y ) - ( P t + u P x + v P y ) = 0, (3. 2d) 


p ( Fj , + u Fj x + v F j y ) = 0, j = l, 2, ... N , 

£ (Fj / Wj ) 

p =P r ~ir 

I (Fj.-IWj) 

. > =1 


(3.2e) 


(3.2f) 


where p is the density, u the velocity component in the x direction, v the velocity com- 
ponent in the y direction, T the temperature, P the pressure, F j the mass fraction of 
species j with molecular weight Wj , and H the mixture enthalpy defined by 

H = ZHj Fj / , (3.2g) 

(=1 


where Hj is the enthalpy of species j 

T 

Hj =Hf I H„ + J C pJ (T ) d T. 

'J'O 


(3.2h) 
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Here C p j is the dimensionless specific heat at constant pressure for each species and the 
subscript °° indicates those quantities that are to be measured in the free stream at 
y = + oo. Finally, the constant C j appearing above is given as 


Ci = 



(3.3) 


The equations were nondimensionalized by selecting the free stream values p^, u^, 
Poo, Ho,, at y = 4* oo for the density, velocity components, temperature, pressure, and 
enthalpy, respectively. Lengths are referred to L, some characteristic length scale of the 
mean flow, while the time scale is the ratio of the length and speed scales. In writing 
down these equations the only assumption made is that all body forces are neglected. 

We assume that the mean flow is given and satisfies the boundary conditions, 

u = p = T = P = l, at y = + ®o , (3.4a) 


u= Pt/, p = P p , T = p Ti P=l aty=- oo, (3.4b) 


where P^ , P p , and Pj are the ratios of the velocity, density, and temperature, respectively, 
in the stream at y = - to that of the stream at y = + «>. Note that P p is related to p^ 
through the identity 


1 = Pp Pt 



(Fj.— 'Wj) 


(Fj'oo/Wj) 


(3.5) 


The stability problem can be formulated independently of the detailed form of the 
mean flow profiles. The mean flow field is perturbed by introducing two dimensional 
wave disturbances of the form 

[ u, v, P, p, T, H, Fj ] (x, y, t) = [U,0, l,p,T,H,Fj](y) 

+ e [/, a <i>, n, r, h,fj ] (y) (3.6) 

where e < < 1 , a is the complex wavenumber and £b is the frequency which is taken to be 
real because we are only treating the spatial stability problem. Substituting (3.6) into the 
inviscid compressible equations (3.2) yields ordinary differential equations for the pertur- 
bation amplitudes. To be compatible with the mean flow equations, the perturbation 
amplitude equations are rewritten in terms of the similarity variable r\, defined by 

r =b dy ’ "-Ts- 


(3.7) 
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In terms of the similarity variable T|, the linearized equations governing the perturba- 
tion amplitudes are, 

i (U - c )r+pp'0 + p(// + p <(/ ) = 0 , (3.8a) 

p [ i ( U - c ) f + p U 0 ] = -/ Cj n , (3.8b) 

i ( U - c ) <(> = - C ! n / a 2 , (3.8c) 

p [ / ( U - c )h+pH$] = i ( U - c )n, (3.8d) 

i ( C/ - c ) fj + p Fj <|) = 0 , j = 1,2, ... N , (3.8e) 

n = £ [pT fj +T Fj r + P F j %], (3.8f) 

j = i "7 . 

where 

^ = £ [Hjfj+Fj C pJ ^]/H co , (3.8g) 

y=i 

C 2 _1 = £ Fj^IWj , (3.8h) 

7=1 

and 

oo = co / 2 V*" , & = a / 2 Vx" , <> = 2 Vjc~ <b . (3.9) 

Here, c is the complex wave speed 

c = co / a . (3.10) 

Since a is complex, the real part of a is the wave number in the x direction, while the 
imaginary part of a indicates whether the disturbance is amplified, neutral, or damped 
depending on whether a, is negative, zero, or positive, assuming positive group velocity. 
The phase speed, c ph , is given by co / a r . If a, is zero, c = c N is the phase speed of a 
neutral mode. Finally, primes indicate differentiation with respect to the similarity vari- 
able rj. 

From this set it straightforward to derive a single equation governing fl, given by 
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(U -C ) 2 


77-77- -7 1 - n = 0 , ( 3 . 11 ) 

( ~'2 U \ U 2 


where 


G i(T|> = Z [Fj/Wj] , 


(3.12a) 


G 2 0l)= I C„ >7 - ///J , 
j = 1 


(3.12b) 


G 3 0 l) = I F////J . 

;'=i 


(3.12c) 


For a single species gas and constant specific heat, (3.11) reduces to that of Jackson and 
Grosch (1989a). 

The boundary conditions for n are obtained by considering the limiting form of 
equation (3.11) as rj — > ± 00 . The solution to (3.11) is of the form 


n exp ( ± Q ± T| ) , 


(3.13) 


where 


Q+ = a 2 1 - 


(1-c) 

C, 


(3.14a) 


Q 2 = 


— Pt - 

Pp 2 Pr 


(Pt/ -c? 


Ci G 


2 KJ \,-< 


. (3.14b) 


We define c+ to be the values of the phase speed for which Q+ vanishes. Thus, 


c + = 1 


(3.15) 


- Pt/ + VPt" 


(3.16) 


c 2 g x . 
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Note that c + is the phase speed of a sonic disturbance at T) = + 00 and c_ is the phase 
speed of a sonic disturbance at T| = — °°. If we define the mixture quantities 

C p = £ Fj C pJ , (3.17a) 

7=1 

-1 

(3.17b) 

and use the identity for the sound speed 

a 2 = (y- 1) C p T = yP / p , (3.17c) 

then c+ can be rewritten in terms of the mixture Mach number at r\ = +«= 



’ N 


N 

Y = 

lyjFji Wj 


S Fj IV/; 


,;=i 


.M 


as follows 


where 


At 


M =u oa ! a„. 


C+ ~ 1 M 


c_ - + 


M ’ 


Pv = Y-~ / Y~ 


M = M* s 


1 + VP7 \ Pp 


1-Pt/ ’ 


(3.18) 

(3.19) 

(3.20) 

( 3 . 21 ) 

(3.22) 


c ± are equal. Thus, from (3.1), 


M _ M d-Pt/) 
M* 1 + VPy ~i Pp 


(3.23) 


The nature of the disturbances and the appropriate boundary conditions can now be 
illustrated by reference to Figure 3, where we plot c+ versus M c . In what follows we 
assume that a 2 r > a 2 ,-. These curves divide the c r -M c plane into four regions, where c r 
is the real part of c. If a disturbance exists with a M c and c r in region 1, then Q. 2 + and 
Q 2 _ are both positive, and the disturbance is subsonic at both boundaries. In region 3, both 
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Q 2 + and n 2 _ are negative and hence the disturbance is supersonic at both boundaries. In 
region 2, G 2 + is positive and Q 2 _ is negative, and the disturbance is subsonic at +<» and 
supersonic at ->», and we classify it as a fast mode. Finally, in region 4, £2 2 + is negative 
and £l 2 _ is positive so the disturbance is supersonic at -h» and subsonic at -<», and we 
classify it as a slow mode. 

To complete the stability problem, one can now see that the appropriate boundary 
condition for either damped or outgoing waves in the moving and stationary streams are, 
respectively, 


n -► e~ a '\ 

if C r > C + , 


if c r < c + . 

(3.24a) 

n -> e a A 

if C r < C_, 


if c r > c_. 

(3.24b) 


3.2. Results. The effect of the scaling parameter M c can now be analyzed. In our 
previously reported calculations (Jackson and Grosch, 1989b) we considered only a single 
species gas with Py = 0 and p y = 1. For this case the convective Mach number (3.23) 
reduces to 


M, 


M 

i + VPT 


(3.25) 


In order to present the variation of the maximum growth rates with M c we normalize them 
by defining the normalized growth rate by 


(I 

(~ a i )max (Pr-0) 


(3.26) 


Figure 4 is a plot of R versus M c for the three thermodynamic models and Pj of 0.5, 1, 
and 2. This data is taken from Jackson and Grosch (1989b). It can be seen that with 
these scalings the data collapses onto essentially a single curve for M c < 1, and a narrow 
band for M c > 1. One should also note that the second unstable modes appear around 
M c = 1. This curve is similar to that obtained by Ragab and Wu (1988) who use 
Bogdanoff’s heuristic definition of the convective Mach number. However, their graph is 
incomplete and absent of the second mode. 

It is instructive at this point to compare our convective Mach number (3.23) to that 
derived by other authors using hueristic arguments. There are basically three definitions 
used. These are, in our notaion, that of Bogdanoff (1983) 

M(l-P„) 

C " (1 + 1/VPJ) Py /4 ’ 


(3.27) 
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Papamoschou and Roshko (1986, 1988) 

M c = 


M Oj-M 

i + WP^ ’ 


and Zhuang, Kubota, and Dimotakis (1988) 

M c = M (1 -c ph ). 


(3.28) 


(3.29) 


where c ph is the phase speed of the most unstable mode from linear stability theory. For 
a single species gas Py = 1 and our definition of the convective Mach number and that of 
Bogdanoff and Papamoschou and Roshko are equal. For most multi-species gases Py is 
close to one and all three definitions yield values of the convective Mach number that are 
about the same. Table 1 contains a comparison of the values of the convective Mach 
number using the definitions given in the equations (3.23), (3.27), and (3.28), and the data 
of Papamoschou and Roshko (1988). The definition of (3.29) requires knowledge of the 
values of the phase speed c ph of the most unstable wave before one can calculate M c . 

Further understanding of the role of the convective Mach number can be gained by 
noting M c = 1 implies that M = M*. This is the largest value of the Mach number for 
which any subsonic instability waves can exist. For larger values of the Mach number 
there are only supersonic modes which radiate into one or the other stream. The largest 
value of the convective Mach number for which only subsonic modes can exist is given 
by 

M Cm , n ~ min(M + , M_) / M* , (3.30) 


where M* is given by (3.22) and M + is the value of the Mach number for which c + is 
zero, and M_ is the value for which c_ is one. In the intermediate range M c < M c < 1 
both subsonic and supersonic instability waves can exist. 

All previous results are for two dimensional disturbances. If three dimensional dis- 
turbances were allowed, the only effect would be to replace the Mach number M by the 
effective Mach number M cos (0), where 0 is the direction of the propagation of the dis- 
turbance. Thus, an effective convective Mach number can be defined by 

M Ctff =M c cos(Q), (3.31) 

with M c given by (3.23). The introduction of the cas(0) factor to produce an effective 
convective Mach number has been suggested by a number of other authors based on 
heuristic arguments. Figure 5 contains curves of the normalized maximum growth rate 

R _ (~ a i )max (P 

( _a ( )max (Pr,0,0°) 


(3.32) 


1 
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versus the effective convective Mach number for the Tanh model with pj = 1/2 and 
P^ = 0. These curves are based on the numerical calculations of Jackson and Grosch, 
1989a. Results are shown for 0 = 0°, 20°, 40° , 60° . For angles of propagation greater 
than 60° it was found (Jackson and Grosch, 1989a) that the maximum growth rate begins 
to decrease with increasing angles of propagation. The data points shown in this figure 
are taken from the experimential results of Papamoschou and Roshko, 1988. Morkovin 
(1987) has suggested that three dimensional disturbances dominate the flow field once the 
effective convective Mach number exceeds one. This insight is borne out by the results 
shown in this figure. 

The above analysis is strictly valid when is not close to one. If P^ were taken to 
be one, then U = 1 is a solution for the mean velocity profile, but the temperature would 
still vary throughout the flow field (if P 7 - is not equal to one). Any instability present 
would be a purely thermal one and the physics would be different. Thus decreasing M c 
by increasing p y to one is inconsistent with the formulation of the problem. It appears 
that the only ways to decrease the convective Mach number is to increase the ratio P 7 / P p 
or to increase the angle of propagation of the disturbances. 

Finally, the above formulation is also valid for a mixing layer in a channel with a 
zero streamwise pressure gradient. 
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Gas 1 

Gas 2 

M 

Pp 

Pi/ 

K 

(3.23) 

M c 

(3.27) 

M c 

(3.28) 

At 

Ar 

3.4 

0.43 

0.81 

0.26 

0.26 

0.26 

n 2 

n 2 

3.1 

0.54 

0.74 

0.34 

0.34 

0.33 

n 2 

Ar 

2.8 

1.80 

0.75 

0.39 

0.37 

0.39 

n 2 

Ar 

3.2 

1.20 

0.67 

0.53 

0.53 

0.55 

He 

N 2 

1.7 

9.20 

0.52 

0.60 

0.61 

0.64 

At 

Ar 

3.1 

0.24 

0.13 

0.89 

0.89 

0.89 

He 

n 2 

2.6 

5.50 

0.42 

1.03 

1.01 

1.05 

He 

n 2 

3.4 

2.20 

0.29 

1.39 

1.38 

1.44 

He 

At 

3.1 

2.40 

0.04 

1.81 

1.81 

1.81 


Table 1. Comparison of the convective Mach number values using the 
definitions of this paper (3.23), Bogdanoff (3.27), and Papamoschou and 
Roshko (3.28) and the data of Papamoschou and Roshko, 1988. 
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Figure 3. Plots of the sonic speeds c ± versus convective Mach number. 
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Figure 4. Plot of the normalized maximum growth rate as a function of the 
convective Mach number for the various thermodynamic models and values 
of pj- taken from Jackson and Grosch, 1989b. 
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